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Abstract. We study the modification of the multiplicity distributions in MLLA 
due to the presence of a QCD medium. The medium is introduced though a 
multiplicative constant (fmed) in the soft infrared parts of the kernels of QCD 
evolution equations. Using the asymptotic ansatz for quark and gluons mean 
multiplicities (nc) ~ e"'^ and (ng) = r~^e^^ respectively, we study two cases: 
fixed 7 as previously considered in the literature, and fixed a^. We find opposite 
behaviors of the dispersion of the multiplicity distributions with increasing fmed 
in both cases. For fixed 7 the dispersion decreases, while for fixed Os it increases. 

1 Introduction 

The standard explanation for the phenomenon of jet quenching - the suppression of parti- 
cles with large transverse momentum produced in nucleus-nucleus collisions compared with 
the expectations from an incoherent superposition of nucleon-nucleon ones - observed at the 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory [1], is radiative 
energy loss (see the reviews (2j). This explanation implies a modification of the standard QCD 
radiation pattern (see e.g. [3] for a recent proposal). 

In standard perturbative QCD, the process of gluon radiation suffers from several types 
of singularities at each order in a^, see Fig. [1] leading to single and double logarithms, which 
have to be resummed. More specifically, there can be two different singularities when a gluon 
is emitted: collinear singularities when the gluon emission angle is very small {6 0) leading 
to ln0; and infrared divergences when the emitted gluon takes a very small fraction x of the 
energy of the parent, leading to ln(l/a;). These singularities give problems of convergence, thus 
the need of resummations. There are different resummation schemes, namely: 

— Leading Logarithmic Approximation (LLA): it resums the single logs [a^ In with 
kr the transverse momentum of the emitted gluon with respect to the parent parton. 
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Fig. 1. Kinematical variables in the process of gluon radiation from a parent parton. 
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— Double Logarithmic Approximation (DLA) : resummation of coUinear and infrared singular- 
ities [asln(l/a;)ln6l]" ~ 

— Single Logarithms (SL): takes into account the emission of hard coUinear gluons [6 — > 0), 
[as ln0]" -O(V^). 

— Modified Leading Logarithmic Approximation (MLLA): SL correction to DLA, 
[a,ln(l/a;)ln6l + a,ln6l]" - [0{l) + {0{^s)]- 

The standard radiation pattern in QCD [4] has been extensively tested by measurements 
in high energy e"'"e~ and pp{pp) collisions. While the DLA resummation scheme is known to 
overestimate the cascading process, as it ignores the energy-momentum conservation since the 
radiating particle does not loose any energy after soft gluon radiation, it has been found a 
very good agreement between the data from OPAL and TASSO and the MLLA approximation, 
while TeVatron data require further refinements [5J . 

When a medium is created in a heavy ion collision, the radiating particle will travel through 
it. A distortion in the radiation pattern will appear. In this work, we focus on the multiplicity 
distributions within a jet. We will use a simple medium implementation within MLLA as 
proposed in [6], to describe the QCD radiation pattern in medium. We go beyond previous 
studies [7j by considering different possibilities to fix the parameters in the calculation, which 
produces qualitivavely different results. 



2 Calculation and results 

The computation of multiplicity distributions in QCD has been discussed many times, see e.g. 
[4l8l9j . Here we simply indicate the main formulae: The multiplicity distributions at a scale 
(rapidity) y — hi(p9 /Qq) — ln(2(5/(5o)i with Qq an infrared cut-off, is the probability of finding 
n objects in the final state radiated within an angle 9 from a parton with momentum p: 

P^y) = = ^ . (1) 

The generating functional is defined in the fohowing way: 

oo 

G{y,z) = Y,P^{y){l + zr. (2) 



The normalized factorial moments of range q are defined as 



E„Pnn(n-l)...(n-g+l) 



2 = 



1 



(3) 



so the generating functional can be expressed as a function of those moments: 



g=0 

The normalized factorial moment of range one is equal to 1. The normalized factorial moment 
of range two, 

{ny 

is related to the normalized dispersion of the distribution. 
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in the following way: 



(n) 



The master evolution equation for splitting of a parton A into partons B and C, with 
momentum fractions x and (1 — x) respectively, is [4j 



'^y - B.c 



dx^^^'i'- [GB{y + Inx, z)Gc{y + ln(l - x), z) - Ga(2/, z)]. (8) 







The different possibilities for splitting [their respective probabilities] are: g g{x)g{\ — x) 
['Pf'ix)], g ^ q{x)q{l-x) [<?g^(x)], q ^ q{x)g{l - x) [<?§^(x)] and q ^ g{x)q{l - x) 

[<?q'^(x) = ^q'~^(1 — x)]. The evolution equations for the generating functional of gluon and 
quark jets in MLLA read 



1 



G'G{y,z)= dxKS{xh^o[GG{y + lnx,z)GG{y + \n{l~x),z)~GG{y,z)] (9) 







+ nf / dxK2{x)j^[GQ{y + lnx, z)GQ{y + \n{l - x),z) ~ GG{y,z)], 







1 



G'Q{y,z)^ / dxK^{xh^o[GG{y + lnx,z)GQ{y + ln{l-x),z)-GQ{y,z)], (10) 







where the prime denotes the derivative with respect to y, and 70 — 2Ncas/TT, with Nc = 3 
the number of colors and nf{= 3) the number of active flavors. The negative terms in both 
equations take into account the probability of no radiation and yield the corresponding integral 
finite without the need of any regulation of the splitting kernels. 

In the absence of any medium, the splitting kernels Kl are the standard DGLAP kernels 
e.g. used for jets in e^e~ annihilation. We introduce the medium in the splitting kernels by 
multiplying the infrared-divergent parts by a constant fmed as done in [6]: 



K${x) = i±i™££ - (1 - a;)[2 - x{l ~ x) 



(11) 



^G( ^ Gf 



Nr. 



1 ~1~ fmed 



where Cf = {N^ ~ 1)/ (2iVc) , and represents the splitting of a particle A with momentum p 
into a particle B with momentum fraction xp and another particle C with momentum fraction 
(1 — x)p. For fmed = wc obviously recover the standard DGLAP kernels. The constant fmed 
grows with increasing medium length and transport coefficient, and decreasing parent-parton 
energy, but the quantitative relation has not been established yet. 

Now we proceed by substituting ^ into ^ and ^TU\i . and by collecting terms with equal 
powers of the auxiliary variable z in both sides of the equations. In this way we will obtain 
equations for the mean multiplicities at order z, and for the corresponding normalized factorial 
moments of range two at order z'^. 
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2.1 First order equations: 0{z) 
At lowest order in z wc get: 

{nG{y)y= I dx^lK^{x)[{nG{y + \^x)) + {nG{y + H^~x)))-{nG{y))] (12) 



+ Uf / dxY^K'^{x)[{nQ{y + Ina;)) + {nqiy + ln(l - x))) - {naiy))], 
Jo 

My))' = / dxj^K^{x)[{nG{y + Inx)) + {ngiy + ln(l - x))) - {ngiy))], (13) 



in which Ui has the meaning of the parton mutiphcity within a jet generated by a parton of 
type i = G,Q. 

Using the asymptotic ansatz for the mean multiplicities [4|8j : 



(ug) = {ng) = ie^^ (14) 

which is valid at very high energies, we get from and 

7 = / dxj^ \kS{x)[x-' + (1 ~ a:)T - 1] + '^K^{x) [x'< + (1 - x)'' - r]l , (15) 
Jo ^ ^ -I 

7=/ dx72xg(a;) [ra;'' + (1 -x)"^ - 1] (16) 



respectively. By doing that, we have obtained two algebraic equations instead the original 
integro-differential ones. The only variable is the medium characterized by fmed, and we have 
a system of two equations with three unknown parameters. Fixing one of them we can solve 
the system analytically. We study two cases: 

— We fix 7 as done in [7] and we obtain the dependence of as and r with fmed- 

— We fix as and study the variation of 7 and r with fmed- 

The results are shown in Figs. [2] and [3l We fix the values of 7 and as in order to agree, for 
fmed = 0, with the results in [7]. 

When 7 is set to be medium independent, which with the ansatz (jl4[) means that the parton 
mean multiplicity does not change when the jet travels through a medium of increasing density, 
a decreasing behavior of as is found, with — > when fmed 00. In Fig. [^bottom we 
indeed see a very drastic decrease of ag in a range of fmed from to 4. As fmed is expected to 
be proportional to g cx T"^ 6 10 , we find a large decrease in the coupling constant for a small 
increase in temperature. The ratio of gluon- and quark-jet mean multiplicities, r, increases with 
the medium (Fig. [3]), with the large fmed limit 



r — > 7 



[lE + f (7)] ( 1 - ^) + ^ 



(17) 



meaning that even in the medium gluon jets are still more active than quark jets in producing 
secondaries. 

In the second case wc analyze, we fix as and calculate the dependence of the anomalous 
dimension 7 with the medium. As shown in Fig. [l^op, it increases. The large fmed limit is 

7, ^ 1. (18) 

Thus, the mean multiplicity of partons increases with the presence of a medium. When we 
calculate r, we find a decreasing behavior (Fig. [3]) with an asymptotic value of 1 when fmed — *■ 
cxD. So the quark and the gluon jets tend to produce secondary partons in the same amount 
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when traversing a medium. This result suggests that for large medium effects, multiplicity is 
dominated by the gluon piece. It agrees with the behavior observed in of r decreasing with 
increasing 7. But it does not coincide with the values r = Nc/Cp = 9/4, 7 = \/7o(l + fmed) 
obtained when taking just the infrared-divergent parts of the splitting kernels. It is actually 
driven by the non-infrared divergent pieces which cancel for x — *■ in (jlSp and (jl6p and, at 
first glance, it looks dubious. It might be due to the ansatz We leave this point for future 
investigations. 

After obtaining the parameters by solving the equations 0{z), we substitute them in the 
equations C(z^). 



2.2 Second order equations: 0{z^) 



Collecting terms with powers in both sides of equations ^ and pUj) we get: 



^{{nQ{y)ry ^ dxilK<^{x) (19) 
X [i(7iG(y + \nx)fF2 + \{nQ{y + ln(l - x))fQ2 
+ {naiy + \^x)){nQ{y + ln(l - x))) - ]^{nQ{y)fQ2l 
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r=<nQ>/<nQ> 



fixed 0^=0.1156 
fixed -1^0.3746 
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Fig. 3. Ratio between gluon and quark mean multiplicities versus fmed, for both cases of fixed as = 
0.fl56 (red solid line) and fixed 7 = 0.3746 \7\ (blue dotted line). 

^[{na{y))^)' ^ j\x^lK'i{x)\ (20) 

X [{{nG{y + \nx)f + {naiy + ln(l - x))f) F2 

+ (noiy + liix)){nG{y + ln(l - x))) - i(nG(j/))'i^2] 
rl 1 

+ nf dxj^K§{x)-[{{nQ{y + \nx)f + {ngiy + \n{l - x)))^) Q2 

+ {ngiy + \nx)){nQ{y + ln{l - x))) - ^{nQ{y))^F2] 

for the second order normalized factorial moments for quarks {Q2) and gluons {F2) respectively. 
Note that in MLLA the second order normalized factorial moments are y-independent. 
Introducing our ansatz (jl4[) for the mean multiplicities, these equations read 

2Q27= / dx-f^Kg{x){r^x^^F2 + {l-xf^Q2 + 2rx^{l-xy ~Q2) , (21) 



27^2= / dx-f^Kg{x){x^''F2 + {l-xf''F2+x^l^xy ^F2) 
Jo 

+ ^ [ dx^lK^ (x2tQ2 + (1 - xf^Q2 + x^{l- xy - r^F2) 



(22) 



As we know the analytic expression of as{fmed) in the case of fixed 7 or ^{fmed) for the 
case of fixed as, and r(fmed), the only unknowns are Q2 and i^2- Solving these system of two 
algebraic equations we obtain F2{fmed) and Q2{fmed)- We show the results in Fig.lH where we 
can see a decrease in the second order moments for both quark and gluon jets. The decreasing 
behavior is found for the two different cases that we analyze. 
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Fig. 4. Top: second order normalized factorial moment Q2 for quarks versus fmed- Bottom: second 
order normalized factorial moment F2 for gluons versus fmed- The line convention is the same as in 
Fig.El 



In order to see what really happens to the multiphcity distribution with the medium, we 
calculate the dispersion by using the relation given by ([7]): 

Dl = {uQ^Dl^^,,^ = ;^e2^'(Q2 - 1 + re-^), (23) 
Dl = {nG?Dl„ = e'^y{F, 1 + g-^^). (24) 

While, as indicated above, in MLLA the second order normalized factorial moments are y- 
independent, this is not the case for the dispersion. On the other hand, the ansatz (fT4|) is only 
valid at large rapidities. Therefore, for the purpose of illustration we choose a large rapidity 
y = 10.72 which roughly corresponds to the total rapidity interval at top RHIC energy. 

In Fig. [5] we show the results obtained for the dispersion of the multiplicity distributions of 
quark and gluon jets. Now the behaviors in the two cases that we analyze are very different: 
For fixed 7 we obtain a decreasing dispersion with increasing fmed- For fixed as an increase 
in the dispersion of the multiplicity distributions is obtained. The former result is difficult to 
understand, while the latter makes more sense within the logic that the multiplicity distributions 
should get wider in medium than in vacuum. 

On the other hand, the normalized dispersions can be directly read from ([23|) and p4|) . 
For the large rapidity we are examining, the results are similar to those for the second order 
moments in Fig. [4] shifting the vertical scale down by one unit - the correction due to the last, 
rapidity-dependent term in the r.h.s. of the equations is smaller than +0.04. 



Will be inserted by the editor 



Dispersion (y= 10.72) 



^ 100000 




0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 4 



'med 



med 



Fig. 5. Dispersion for gluons (plots on the left) and quarks (plots on the right) versus fmed, for both 
cases of fixed as — 0.1156 (red solid line, upper plots) and fixed 7 = 0.3746 (blue dotted line, lower 
plots), for y = 10.72. 



3 Conclusions 



Analytical estimates of the medium modification of parton branching arc most needed, both 
by their own interest and as a tool for understanding results which come from Monte Carlo 
simulators [11]. These medium modifications will hopefully be tested both at RHIC (see the 
first jet results in [12]) and at the Large Hadron Collider [13] in the near future. 

In this work, we have considered the modification of the multiplicity distributions in MLLA 
due to the presence of a QCD medium. The medium has been introduced though a multiplicative 
constant (/med) in the soft infrared parts of the kernels of QCD evolution equations. Using 
the asymptotic ansatz (1141) for quark and gluons mean multiplicities, we have studied two 
cases: fixed 7 as previously considered in [7], and fixed as- After analyzing several features 
of both cases, we have discussed the mutiplicity distributions. We find opposite behaviors for 
the dispersion of the multiplicity distributions with increasing fmed- For fixed 7 the dispersion 
decreases, while for fixed as it increases. As future developments, we plan to clarify several 
points in our analysis, both relaxing the ansatz (|14p and including medium effects in the splitting 
functions with a more refined treatment than the use of fmed- 
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